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The theory of hole superconductivity predicts that in superconductors the charged superfluid is 
about a million times more rigid than the normal electron fluid. We point out that this physics 
should give rise to large changes in the bulk and surface plasmon dispersion relations of metals 
entering the superconducting state, that have not yet been experimentally detected and would be 
in stark contradiction with the expected behavior within conventional BCS-London theory. We also 
propose that this explains the puzzling experimental observations of Avramenko et al[T] on electron 
sound propagation in superconductors and the puzzling experiments of W. de Heer et al[2] detecting 
large electric dipole moments in small metal clusters, as well as the Tao effect[3] on aggregation of 
superconducting microparticles in an electric field. Associated with the enhanced charge rigidity is 
a large increase in the electric screening length of superconductors at low temperatures that has 
not yet been experimentally detected. The physical origin of the enhanced charge rigidity and its 
relation to other aspects of the theory of hole superconductivity is discussed. 

PACS numbers: 



I. INTRODUCTION 

A normal metal screens electrostatic fields over dis- 
tances of order of the interelectronic distance, or hp 1 , 
with hp the Fermi wavcvector, quantitatively of order 
A~ l for normal metallic densities. Conventional BCS- 
London theory predicts that the response of supercon- 
ductors to a static electric field is essentially the same as 
that of normal metals [4H9] . Instead, the theory of hole 
superconductivity[in] predicts that superconductors can 
only screen electrostatic fields over much larger distances, 
in their ground state over distances of order Al, the Lon- 
don penetration depth, quantitatively of order hundreds 

of Ami- 

Within the theory of hole superconductivity the inabil- 
ity of superconductors to screen over smaller distances 
originates in the fact that superconducting electrons re- 
side in highly overlapping orbits of radius 2Al[12], and 
this precludes the possibility of charge fluctuations of 
shorter wavelengths that would destroy the ability of the 
superconducting electrons to maintain phase coherence 
as they traverse these large orbits. That superfluid elec- 
trons reside in orbits of radius 2Al follows from the fact 
that according to this theory this is the only consistent 
way to explain dynamically the Meissner effect exhib- 
ited by all superconductors^^. The inability of super- 
fluid electrons to screen leads, via the compressibility sum 
rule, to the prediction that superfluid electrons are highly 
incompressible unlike normal metal electrons. 

This enhanced rigidity of the superconducting fluid im- 
plies that the longitudinal plasmon dispersion relation 
will be much steeper in the superconducting than in the 
normal state[H]. In contrast, BCS theory predicts no 
change in the plasmon dispersion relation 5 8 . The bulk 
plasmon dispersion relation can be measured by EELS 
(electron energy loss spectroscopy) [HI [TS] and by inelas- 
tic X-ray scattering [17] as well as optically in transmis- 
sion experiments through thin films [18] [19]. Also the sur- 
face plasmon dispersion relation, which can be measured 



in EELS [2(3 or optical experiments [2T]. should change in 
the superconducting state. To our knowledge these ex- 
periments have not been yet done on superconductors. 
We discuss here what we expect the observations will 
show, in stark contrast with what would be expected 
within BCS-London theory. 

We furthermore discuss three experiments that have 
been performed in recent years that provide strong evi- 
dence in favor of the enhanced rigidity of superfluid elec- 
trons predicted by our theory: (i) sound propagation by 
electrons (Avramenko effect) [1], (ii) electric dipole mo- 
ments of small metal clusters (de Heer effect) [2] , and (iii) 
aggregation of superconducting microparticles in large 
electric fields (Tao effect) [3]. 

The larger electric screening length of superconduc- 
tors should be directly detectable experimentally. So far 
the only experimental indication of this appears to be a 
report by Jenks and Testardi[22j that measured an in- 
creased penetration of electric field in YBCO films be- 
low T c . We discuss the expected behavior of the electric 
screening length below T c within our theory. 



II. ELECTRODYNAMIC EQUATIONS FOR 
SUPERCONDUCTOR 



Within the theory of hole superconductivity the first 
London equation for the time derivative of the supercur- 
rent is modified to read|14j 



dt 



(E + V</>) 



(1) 



with <f> the electric potential (the V</> -term is absent in 
the conventional London equations). The magnetic vec- 
tor potential A in the second London equation 
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obeys the Lorenz gauge V • A = —(l/c){d(j>/dt) rather 
than the London gauge V • A = OpJ]. Note that Eq. 
(1) follows from Eq. (2) on using Faraday's law. The 
charge density in the superconductor p(r, t) satisfies the 
equation [H] 



dt 2 



(p - Po) 



c 2 V 2 p 



(3) 



and the electric potential cj>(r, t) satisfies the same equa- 
tion 



d 2 (j) 
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where po is a uniform positive charge density and 4>o(r) 
is the resulting electrostatic potential (V 2 </>o = —Anpo). 
The London penetration depth Az, is given by the usual 
formCS] 
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with ujp the plasma frequency. The parameter po is de- 
termined by the condition that the internal electric field 
that develops in the interior of the superconductor due 
to expulsion of negative charge to the surface [TT] should 
reach its maximum value 124 



E„ 



he 
4eAl 



(6) 



within a London penetration depth from the surface, 
pointing outward perpendicular to the surface. 



III. DIELECTRIC FUNCTION AND 
COMPRESSIBILITY 

It follows from the electrodynamic equations discussed 
in the previous section that the longitudinal dielectric 
function for the superfluid within our theory is given 

byna 



e s (k,uj) = 1 - 



c 2 k 2 



(7) 



Eq. (7) is of the generic form of a hydrodynamic longi- 
tudinal dielectric function for the electron fluid 25l [Ml 
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that yields for the static dielectric constant 
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with k the electronic compressibility. The second equal- 
ity follows from the compressibility sum rule|27j. so that 
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with 
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For the free electron gas, the zero temperature com- 
pressibility is 



2nep nm P v\ 
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with tp = m e vp/2 the Fermi energy and vp the Fermi 
velocity, yielding 



(13) 



Instead, for the superconductor we have from Eqs. (7) 
and (8) 



fi 2 = c 2 



(14) 



(15) 



so that the superconducting electron fluid is enormously 
more rigid than the normal metal electron fluid, since 
c >> vp. We expect this enhanced rigidity to show up 
in experiments where electron density oscillations are in- 
duced that are not accompanied by motion of the ions so 
that an electric potential builds up in the interior of the 
superconductor. 

Eq. (9) yields the static longitudinal dielectric func- 
tions for the superfluid electrons and the normal metal 
electrons respectively 
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e s (M) = l + Y2p 



e„(fc,0) = 1 + 
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with Ax, the London penetration depth given by Eq. (5) 
and Xtf the Thomas Fermi screening length given by 
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with ciq the Bohr radius. 



IV. BULK PLASMONS 

The bulk plasmon dispersion relation follows from set- 
ting the longitudinal dielectric function to zero. Eq. (8) 
yields 



(19) 
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hence we predict for superconductors at zero temperature 
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For the normal metal, the plasmon dispersion relation 
obtained from the longitudinal dielectric function calcu- 
lated in the random phase approximation (Linhardt di- 
electric function) yields Eq. (19) with [2"8] 



(21) 



which is slightly different from Eq. (13), valid in the 
low frequency limit 29J. Thus, to reproduce the Linhardt 
dielectric function the hydrodynamic form Eq. (8) has 
to include a variation of (3 2 from low to high frequencies. 

Within a two-fluid description of a metal below the su- 
perconducting transition temperature the electronic pres- 
sure results from the sum of the superfluid and the nor- 
mal fluid pressures, hence 
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and the parameter /3 2 is 



1 



1 , 1 



"(- + -) = 



1 



n n 3 
n 5 



L s 2 



(22) 



(23) 



where n n and n s are the densities of normal and super- 
fluid electrons and we have used the high frequency value 
of /3 Eq. (21) for the normal electron contribution. 

Thus, Eqs. (19) and (23) show that as the tem- 
perature is lowered below T c a sharp increase in the 
slope of the bulk plasmon dispersion relation uJf. ver- 
sus k 2 should be seen. In a two-fluid model one has 
n s = n(l — t 4 ),n n = nt 4 , with t = T/T c , which is also 
approximately the behavior predicted by BCS theory [53], 
so we expect for the bulk plasmon dispersion relation 
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V. SURFACE PLASMONS 

Surface plasmons (also called surface plasmon polari- 
tons) arc longitudinal charge oscillations coupled to an 
electromagnetic wave with both longitudinal and trans- 
verse field components propagating along the surface of 
a metal, excited by either fast electrons or electromag- 
netic radiation. Crowell and Ritchie [3D] and Fuchs and 
Kliewer[5S] derived the following dispersion relation[3"U] 

2r/7 2 2 i 2 2\l/2 , /; 2 2 2\l/2l 

UJ k [{k C + UJ p - UJ k ) I +{k C - L)%) I ] X 

(fc 2 /3 2 +^-^)V 2 = 

^[k 2 ^ + (k 2 c 2 - c^) 1 / W + co 2 p ^)V2 K 25) 
In the limit j3 — > the solution is 
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FIG. 1: Surface (solid lines) and bulk (dashed lines) plasmon 
dispersion relations for different values of j3/c (numbers next 
to the lines). /3 increases as the temperature is lowered (Eq. 
29)). The values of P/c of 0,0.1,0.2 and 0.5 correspond to 
values of T/T c of f , 0.998, 0.990 and 0.93f respectively. 



and UJk — > i^p/V^ for large k. For any j3 ^ 0, the surface 
plasmon dispersion relation for large k is, from Eq. (25) 



= uj 2 p + P 2 k 2 



(27) 



so it approaches the bulk plasmon dispersion relation Eq. 
(19). For small wavevectors the surface plasmon disper- 
sion relation is 
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so it increasingly deviates from the transverse dispersion 
relation uj^. = ck for smaller /3 and larger k. 

Figure 1 shows examples of the surface and bulk plas- 
mon dispersion relations for various /?. As function of 
temperature, /3 increases very rapidly as T is lowered be- 
low T c according to Eq. (23). For any typical value of 
vf (of order 1% of the speed of light) the vf term in Eq. 
(23) can be ignored, so that 



P 2 (T) = (l-t 4 )c 2 



(29) 



The values of (3 of 0,0.1,0.2 and 0.5 shown in Fig. 1 
correspond to values of T/T c of f, 0.998, 0.990 and 0.931 
respectively. Consequently we expect rapid changes in 
the observed surface and bulk plasmon frequencies as the 
temperature is lowered below T c , in contrast to conven- 
tional BCS-London theory that predicts no change 8J. 

However, for surface plasmons the interpretation of ex- 
periments could be more complicated because it appears 
that in experiments performed in the normal state the 
induced charge fluctuations can spill out of the surface, 
drastically modifying the dispersion rclation|3f 33 , an 
effect which is not taken into account by Eq. (25). We 
expect this spill-out effect to be even more pronounced in 
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the superconducting state because of the enhanced rigid- 
ity and because the superconductor has an enhanced ten- 
dency to spill out electrons within our theory [54]. 



VI. PLASMON EXPERIMENTS 

Measurement of the angular dependence of scattered 
electrons in electron energy loss experiments (EELS) pro- 
vides information on the bulk plasmon dispersion rela- 
tion. A large number of such studies has been performed 
on many different normal metals since the 1950's[T5]. 
Generally these studies are done at room temperature, 
although there have also been EELS studies of the effect 
of temperature on the plasma frequency down to liquid 
helium temperatures for ^4i|35| and Pfc|36|. However to 
our knowledge there has not been a single EELS study of 
plasmons in a superconducting metal that would look at 
possible changes in the plasmon dispersion relation below 
the critical temperature (except for ref. |37| for a high T c 
cuprate that did not detect any change presumably due 
to experimental accuracy limitations). This is very sur- 
prising and we hope such experiments will be done in the 
near future. As discussed in Sect. IV we expect a very 
rapid increase in the plasmon energy for fixed wavevector 
as the system is cooled below T c , 

Longitudinal bulk plasmons can also be excited opti- 
cally with obliquely incident p-polarized (parallel to the 
plane of incidence) electromagnetic radiation [151 [H)] and 
the plasmon dispersion relation can be measured. For ex- 
ample, Lindau and Nilsson[19] obtained the bulk plasmon 
ojk for Ag from transmission experiments through thin 
films of different thicknesses of the order lOOA. The ex- 
periment is done at fixed angle of incidence and each film 
thickness gives a small number (2 in this case) of points 
in the dispersion relation. Andcregg et al[38 measured 
the plasmon dispersion relation for K from oscillatory 
structure in the absorption of thin films of varying thick- 
ness from 27 A to 10(M, and were able to extract up to 10 
data points per film. Again, it would seem straightfor- 
ward to do such experiments with superconducting films 
but not a single study has been performed so far to our 
knowledge. We hope such studies will be done in the near 
future. 

Inelastic X-ray scattering (IXS) experiments can also 
provide information on the bulk plasmon dispersion 
relation [391442J . IXS experiments at low temperatures 
have been performed in recent years for example to study 
the physics of liquid and solid fle |43fl4"5] . However no 
attempt has been made to date to study the plasmon 
dispersion relation of metals like e.g. Al [41] in the tem- 
perature range where they become superconducting using 
this technique. 

Surface plasmon experiments on superconductors have 
never been performed to our knowledge, neither EELS 
nor IXS nor optical. With conventional optical meth- 
ods it is complicated to excite surface plasmons because 
of required matching conditions and rough surfaces are 



needed which introduces additional complications [IB] . 
However, recently developed scanning near-field optical 
microscopy techniques [57] provide the possibility to lo- 
cally excite and detect surface plasmons [35] and may al- 
low for detailed studies of the effect of the onset of super- 
conductivity on the surface plasmon dispersion relation. 

Finally, surface plasmons excited in metal nanopar- 
ticles (Mie resonances) [IS] are sensitive to the longitu- 
dinal dielectric function[5nj and thus are likely to show 
interesting changes due to the change in the dielectric re- 
sponse that we predict upon onset of superconductivity. 
Such experiments have never been done with supercon- 
ducting nanoparticles to our knowledge. 

VII. ELECTRON SOUND ANOMALY 

Avramenko and coworkers pQ apply a longitudinal elas- 
tic wave to the surface of a metal and detect an electric 
potential oscillation at the opposite end of the sample. 
They find two types of signals, one propagating at the 
ordinary sound velocity and a much faster one propagat- 
ing at a speed of order the Fermi velocity, which they call 
"electron sound" . When the temperature of the sample 
is lowered below the superconducting transition temper- 
ature the amplitude of the transmitted signals drops pre- 
cipitously. Avramenko et al point out that this behavior 
has no explanation within the conventional theory of su- 
perconductivity. 

According to Avramenko et al the displacement am- 
plitude at the receiving interface for the electron sound 
signal is 

s 

ues ~ u (30) 

Veff 

where s is the sound velocity and v e ff the velocity of 
electron sound propagation which Avramenko et al as- 
sume is the Fermi velocity vp . u a is the amplitude of the 
elastic vibrations at the interface where the signal is gen- 
erated, ues determines the electric potentials measured 
at the receiving interface ips and tpES f° r sound and elec- 
tron sound. Both potentials decrease precipitously as the 
temperature is lowered below the superconducting T c . 

If the superfluid is very rigid compared to the normal 
fluid as predicted by our theory (Eq. (15)) it is natural 
to expect that the amplitude of longitudinal charge os- 
cillations will rapidly decrease as the temperature is low- 
ered below T c and the superfluid concentration increases. 
Following the behavior of the bulk modulus Eq. (23) we 
argue that the electron sound velocity v e ff in Eq. (30) 
below T c can be estimated by 

V n n V v f 

within a two-fluid description. Figure 2 shows the ob- 
tained behavior of the amplitude of the potentials with 
this assumption, compared to the experimental data of 
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FIG. 2: The decay of the sound signal below T c measured by 
Avramenko et al(T] in Ga for two different values of the sound 
amplitude (jagged solid line (high amplitude) and open circles 
(low amplitude) and theoretical prediction (Eqs. (30), (31)) 
for three different values of Vf/c (smooth solid lines). The 
higher line corresponds to the larger vp/c. The dashed line 
gives the behavior of the Fermi function /(A) (from ref. pQ) 
which would be the expected behavior within the conventional 
theory of superconductivity. 



Avramenko et al for Ga[T] (2004) for three different val- 
ues of vf/c. It can be seen that our curves qualitatively 
and semiquantitavely fit the observations for reasonable 
values of vf/c. An even better fit may result from using 
values of the supcrfiuid concentration derived from mea- 
surement of the temperature-dependent London penetra- 
tion depth rather than the two-fluid model temperature 
dependence assumed here. We argue that the compari- 
son shown in Fig. 2 provides strong evidence in favor of 
the greatly enhanced charge rigidity of superconductors 
predicted by our theory. 



VIII. DE HEER CLUSTER DIPOLE MOMENTS 

In a series of papers, W. de Heer and coworkers^ es- 
tablished that small metallic clusters of Nb, V and Ta 
exhibit a large electric dipole moment at low tempera- 
tures (of order several Debye for clusters of up to 100 
atoms). Through a variety of measurements they found 
very strong evidence that the development of the elec- 
tric dipole moment is associated with the onset of su- 
perconductivity. In contrast, similar clusters of a non- 
superconducting metal, Na, showed essentially no elec- 
tric dipole moments|51j. 

Within our theory a superconducting body expels neg- 
ative charge from the interior to the surface and the re- 
sulting charge distribution is rigid. The distribution of 
electronic charge is determined by the geometry of the 
body and can be obtained by numerical solution of the 
electrostatic equations jUj. Initially we had hoped [TT] 



that the inhomogeneous electronic charge distribution 
predicted by our theory would account for the electric 
dipole moments observed by de Heer et al. However our 
calculations show that the electronic charge distribution 
does not exhibit an electric dipole moment even for sam- 
ple shapes without inversion symmetry [53]. 

However, the distribution of ionic charge in a small 
cluster is discrete rather than continuous, and this fact is 
not taken into account in our calculation. Small metal- 
lic clusters have irregular shapes [MJ [55], in general with 
no inversion symmetry. Gcnerically an electric dipole 
moment will be generated by the ionic charges deter- 
mined by the overall shape of the cluster as well as by 
the discrete location of the ions. In the normal state, 
as well as within conventional BCS theory, metallic or 
superfluid electrons are extremely efficient at screening 
electric fields over a length scale \tf, of order lA (eq. 
(18)) and thus will screen any ionic dipole moment. In- 
stead, within our theory the superfluid electrons can only 
screen electrostatic fields over distances of order the Lon- 
don penetration depth (Eqs. (16), (5)), typically of order 
several hundred A, which is much larger than the linear 
dimensions of the de Heer clusters (which have up to 
~ 100 atoms and linear dimensions smaller than lOA). 
Therefore, we argue that the observation of large un- 
screened electric dipole moments in metallic clusters of 
dimensions much smaller than the London penetration 
depth is strong evidence in favor of the large rigidity of 
the superfluid electron charge distribution predicted by 
our dielectric function Eq. (7). 



IX. TAO EFFECT 



In a series of papers[3j, R. Tao and coworkers found 
that superconducting microparticles in a strong electro- 
static field assemble into spherical shapes of macroscopic 
dimensions. We have proposed a detailed explanation of 
this "Tao effect" [SS] , based on the charge expulsion and 
resulting electric fields in the neighborhood of supercon- 
ducting particles of non-spherical shape predicted by our 
theory. 

However, even without considering the details of our 
theory, in a more general context it is clear, as pointed 
out in the experimental papers [3J, that this observation 
is impossible to explain unless electrostatic fields pene- 
trate the superconducting particles a distance consider- 
ably larger than the Thomas Fermi length. This then 
requires that the superconducing charge distribution is 
more rigid than in the normal state where it can screen 
the electric field beyond an A or so of the surface. Thus, 
we argue that the observation of the Tao effect is also 
a strong indicator that the charge distribution in super- 
conductors is more rigid than in the normal state. 
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X. SCREENING OF ELECTROSTATIC FIELDS 

The electrodynamic equations of our theory 
predict fl4]. according to Eq. (16), that the super- 
fluid electrons screen applied static electric fields over a 
distance Al rather than over a Thomas Fermi screening 
length as predicted by BCS theory [8] [57]. In fact, the 
London brothers themselves considered electrodynamic 
equations for superconductors predicting such behavior 
in an early version of their theory 58J. However, shortly 
thereafter H. London performed an experiment [59] 
attempting to detect this effect and didn't find it, after 
which the London brothers discarded that version of the 
theory and adopted the conventional London equations 
to describe superconductors which do not allow electric 
fields in the interior. 

We expect the electric screening length to increase con- 
tinuosly from Xtf ~ l-A to Xt, ~ 100's of A as the tem- 
perature is lowered from T c to 0. In a 2-fluid model 
description the static dielectric constant at finite tem- 
peratures is given by 



e(*,0) = e s (M) + e„(fc,0)-l 

1 1 1 

1 + ( A L (T) 2 + X TF (T) 2 '¥ 

giving the effective electric screening length Xe as 
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(32) 



(33) 



witn n s (t) = n(l — i 4 ), n n (t) = nt A , with t — T/T c . 
Eq. (33) predicts the temperature dependence of the 
electric screening length shown in Fig. 3. Note that only 
at temperatures well below T c does the screening length 
increase substantially. 

In H. London's 1936 experiment [59] he attempted to 
measure changes in the capacitance of a capacitor with 
superconducting electrodes of the metal Hg that would 
result from an increased electric penetration depth. His 
experiment showed no change, from which he concluded 
that the electric screening length doesn't change in su- 
perconductors. However, the lowest temperature reached 
in H. London's experiment was T — 1.8°K, which cor- 
responds to T/T c = 0.43 for Hg (T c = 4.153° K). With 
the sensitivity of his experiment, London could have de- 
tected a change in the capacitance corresponding to the 
screening length increasing above A^ ~ 20A As can be 
seen in Fig. 3, for T/T c ~ 0.4 the screening length would 
only have increased to about 5A, hence substantially less 
than what could have been detected with the sensitivity 
of that experiment. 

In 1993 Jenks and Testardi attempted to measure 
the penetration of a static electric field into epitaxial 
thin films of Y Ba-zCu^Oi-x 22 . They reported a large 
change close to T c , in apparent disagreement with both 
BCS theory and with our expected behavior shown in 
Fig. 3. However, it is not clear that this experiment 
was free of experimental artifacts, since the results also 
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FIG. 3: Screening length for electrostatic field versus temper- 
ature for three values of the zero temperature London pene- 
tration depth (A_l = 790A, 390A and 160A) appropriate for 
Hg, Nb and Al respectively. The inset shows the same data 
on a different scale. Xtf = 1A is assumed. 



showed a large change in penetration depth with tem- 
perature above T c , and variations between different films. 
The experiment has not been repeated, nor are there any 
other published reports of attempts to measure changes 
in the electric screening length below the superconduct- 
ing critical temperature in either high T c or conventional 
materials to our knowledge. 



XI. DISCUSSION 

The need for reformulation of London electrody- 
namics arose in our theory from the prediction that 
negative charge expulsion occurs in the transition to 
superconductivity |60j . which is a consequence of the mi- 
croscopic physics of electron-hole asymmetry [ST] and in- 
compatible with the conventional London equations that 
assume that no electrostatic fields can exist in supercon- 
ductors. Our reformulation renders the theory relativis- 
tically covariant|14|. unlike conventional London elec- 
trodynamics, and allows for a natural and consistent 
extension of the electrodynamics equations to the spin 
sector [24 so as to describe both charge and spin cur- 
rents, which is necessitated by the predicted existence of 
an outward pointing electric field in the interior of super- 
conductors. 

The enhanced charge rigidity and inability to screen 
can be seen also as a natural consequence of several other 
aspects of the theory. For example, superconductivity in 
this theory is driven by kinetic energy rather than poten- 
tial energy lowering [6"2l 155] . Thus, in contrast to the nor- 
mal metal the superconductor is willing to pay a price in 
Coulomb potential energy in order to optimize kinetic en- 
ergy, which naturally results in its inability to effectively 
screen electric fields over short distances, a process which 
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is potential-energy driven in the normal metal. Kinetic 
energy lowering is associated with the fact that in the 
transition to superconductivity electrons 'undress' from 
the electron-ion interaction, expand their wavelength and 
no longer "see' the discrete ionic potential[Sl], hence are 
unable to screen perturbations on interatomic distance 
scales as normal electrons do. 

Associated with the much larger screening length is 
the fact that the compressibility of the superfluid is enor- 
mously reduced compared to the normal metal. This is 
related to the enhanced quantum pressure of the super- 
fluid compared to the normal fluid, which is manifest 
in superconductors in the negative charge expulsion and 
in superfluid 4 He in the fountain effect [65] . It does not 
mean however that the pressure is increased by the same 
factor as the rigidity (bulk modulus). We have for the 
superconductor (Eq. (15)) 

B = — = n—— = m e c 2 n (34) 
K On 

and integrating we obtain 

P = m e c 2 (n — no) (35) 

where no is an integration constant, in contrast to the 
normal metal where 

P = 2 -t F n. (36) 
5 

It is natural to conclude that n—n is of order p_/e (p_ is 
the excess negative charge density near the surface [24]). 
the expelled number density, which is smaller than the su- 
perfluid density n by about the same factor (~ 10 6 ) than 
the energy m e c 2 is larger than the Fermi energy e p |24j . 
Thus the superfluid pressure in the superconductor is of 
the same order of magnitude as the electronic pressure 
in the normal metal, but its rigidity is enormously en- 
hanced. 

How is this compatible with the experimental obser- 
vation that the compressibility of a solid in the super- 
conducting state is essentially the same as in the normal 
state? Clearly in a quasistatic compressibility measure- 
ment the ions and the electrons move together, no charge 
imbalance is generated and the enhanced rigidity does 
not manifest itself. It is only in experiments where the 
electronic density is locally changing relative to the ionic 
density that the much larger rigidity will show up. 

Formally one can write electrodynamics equations for 
the superconductor where the screening length for elec- 
trostatic fields is \l but where no charge expulsion oc- 
curs, as done by the London brothers themselves in the 
early version of their theory 58], as well as by others 
thereafter [66, 67 . Mathematically the formalism is very 
appealing but there is no physics behind it, and per- 
haps for that reason the London brothers were quick to 
discard it soon thereafter when H. London's experiment 
seemed to disprove it[5S]. For us instead, the enhanced 
charge rigidity and the predicted charge expulsion are 




normal Superconducting 
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FIG. 4: Electronic orbits expand from radius kp in the nor- 
mal state (left) to radius 2Al in the superconducting state 
(right). This is the origin of the charge expulsion and the 
charge rigidity over distances of order Al predicted by the 
theory. The black dots denote the instantaneous position of 
the electron, i.e. the "phase", which is random in the nor- 
mal state where the orbits are non-overlapping and highly 
correlated between different overlapping orbits (i.e. phase co- 
herent) in the superconducting state. The orbiting speed is 
v° = h/ '(4m e Ai) in the superconducting state [24]. 



inextricably linked: no enhanced rigidity can take place 
without charge expulsion and no charge expulsion can 
occur without accompanying enhanced rigidity. This is 
because both phenomena are a direct consequence of the 
fact that electronic orbits expand, driven by kinetic en- 
ergy lowering, from microscopic non-overlapping orbits 
of radius kp 1 to orbits of radius 2Al in the transition 
to superconductivity [T2], as shown schematically in Fig. 
4. Orbit expansion implies outward motion of negative 
charge, and the resulting mesoscopic orbits are highly 
overlapping which makes it impossible to create a charge 
fluctuation over a small distance since the extra electrons 
would not have the ability to insert their orbits in the 
mesh of highly correlated interpenetrating orbits that al- 
ready exists. 

It is often said in the context of the conventional the- 
ory of superconductivity that the wavefunction of a su- 
perconductor is "rigid" , a concept first introduced by F. 
London. In the conventional theory, "rigidity" refers only 
to the response to magnetic perturbations. Instead, our 
theory extends the property of rigidity of superconduc- 
tors also to the response to electric perturbations. Rigid- 
ity to both magnetic and electric perturbations originates 
in the overlapping phase-coherent 2\l orbits depicted in 
Fig. 4, which also explains the macroscopic phase coher- 
ence (phase rigidity) of the superconductor: an electron 
orbiting out of phase would collide with other electrons 
in overlapping orbits and pay a high price in Coulomb 
energy. And this also explains why the length Al enters 
symmetrically in our theory for both magnetic and elec- 
tric phenomena |24j: the 2Xl orbits are necessary for the 
Meissner effect to take place [15]. as already suspected 
long ago by Smith[68] and by Slater[69j, and the same 
2A^ orbits determine the electric screening length. The 
wavefunction for the superconducting state has to de- 



scribe superfluid electrons in 2Al orbits, which BCS the- 
ory does not do, if it is to describe the ubiquitous Meiss- 
ner effect and the experimental consequences of enhanced 
charge rigidity discussed in this paper. 

The fact that the superfluid wavefunction is rigid with 
respect to both magnetic and electric perturbations and 
the resulting new electrodynamics follow naturally in a 
relativistic context[SH]. Within Klein-Gordon theory [70] 
describing a relativistic scalar wavefunction ^(r, t), the 
current four-vector J = (J(r,t),icp(r,t)) for the current 
and charge densities is given by 

*(-~V- e A(f,t))y*] (37a) 
i c 

e d 

Vi-ih^-e^t))**} (37b) 

with A the magnetic vector potential and <j) the electric 
potential. In the conventional theory it is said that the 
Meissner effect results from the fact that the wavefunc- 
tion is unaffected by changes in the magnetic vector po- 
tential A. Hence, since J = in the absence of magnetic 
fields, eq. (37a) implies that for any value of A 

2 

_. 77 f>^ _. 

J{r,t) = -^A(r,t) (38) 
m e c 

with n s = ^*^f, giving rise to the Meissner effect. Ex- 
tending the argument it is natural to assume that the 
wavefunction \&(r*, t) is also unaffected by applied electric 
fields and by proximity to the boundaries of the sample. 
If deep in the interior of the superconductor the charge 
density is assumed to be a constant po, with associated 
electric potential 4>o(r) (V 2 0o = — 47rpo), it follows from 
applying Eq. (37b) to a position deep in the interior 



and another arbitrary position r and substracting, that 
at any position r with or without applied electric fields 

2 

77 p 

p(r, t)-p = -^W, t) - Ur)) (39) 

which is the basic equation of our modified electrody- 
namic formalism|14j determining the charge distribution 
and electric potential in superconductors of arbitrary 
shape. 

In summary we have discussed in this paper six differ- 
ent experimental probes of the enhanced charge rigidity 
of superconductors predicted by our theory. Three of 
them (electron sound, de Heer effect, Tao effect) have 
already shown clear evidence for enhanced charge rigid- 
ity. Another two (bulk and surface plasmons) have not 
yet been experimentally tested, and there are a variety 
of different experimental techniques (EELS, IXS, optical 
transmission, optical near-field, nanoparticles) that can 
be used for that purpose. Finally, the predicted increased 
electric screening length below T c has yielded ambiguous 
results so far[22, 59]. For none of the three observations 
that we interpret as arising from the enhanced charge 
rigidity predicted by our theory have alternative plausi- 
ble explanations been proposed, and they all seem to be 
incompatible with conventional London-BCS theory For 
the changes that we predict in the bulk and surface plas- 
mon dispersion relations no other such predictions have 
been made in other theoretical frameworks and they are 
also incompatible with conventional BCS-London theory, 
as is the predicted increase in electric screening length at 
low temperatures. It will be interesting to confront the 
predictions of our theory and of BCS-London theory with 
future experimental results. 
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